ESI-619 
OUT-4102-77 
MZ-TH/98-42 
hep-th/981008T 
13 October 1998 



Renormalization automated by Hopf algebra 

D. J. Broadhurst 1 ) and D. Kreimer 2 ) 

Erwin Schrodinger Institute, A- 1090 Wien, Austria 



Abstract It was recently shown that the renormalization of quantum field theory is 
organized by the Hopf algebra of decorated rooted trees, whose coproduct identifies the 
divergences requiring subtraction and whose antipode achieves this. We automate this 
process in a few lines of recursive symbolic code, which deliver a finite renormalized 
expression for any Feynman diagram. We thus verify a representation of the operator 
product expansion, which generalizes Chen's lemma for iterated integrals. The subset 
of diagrams whose forest structure entails a unique primitive subdivergence provides a 
representation of the Hopf algebra Hr of undecorated rooted trees. Our undecorated 
Hopf algebra program is designed to process the 24,213,878 BPHZ contributions to the 
renormalization of 7,813 diagrams, with up to 12 loops. We consider 10 models, each 
in 9 renormalization schemes. The two simplest models reveal a notable feature of the 
subalgebra of Connes and Moscovici, corresponding to the commutative part of the Hopf 
algebra 1~Lt of the diffeomorphism group: it assigns to Feynman diagrams those weights 
which remove zeta values from the counterterms of the minimal subtraction scheme. We 
devise a fast algorithm for these weights, whose squares are summed with a permutation 
factor, to give rational counterterms. 
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1 Introduction 



Perturbative quantum field theory (pQFT) entails the process of renormalization: the 
iterated subtraction of sub divergences of Feynman diagrams, culminating in the subtrac- 
tion of an overall divergence (or addition of a counterterm) that renders each diagram 
finite and delivers its contribution to the appropriate renormalized Green function. Many 
people helped to develop this process, though it has become customary to record the work 
of Bogoliubov and Parasiuk PL Hepp 0, and Zimmermann || (BPHZ), completed by 
the forest formula in ||, almost 30 years ago. A review is provided by the textbook of 
Collins 0. In addition to the BPHZ formalism one needs analysis, to evaluate regularized 
expressions for the bare diagrams, prior to renormalization. Much progress has been made 
in this direction by combining dimensional regularization || and integration by parts, for 
massless || and massive |7|, § || diagrams with up to 3 loops, together with more ad hoc 
techniques at higher loops. Since 1981 it was possible, in principle, to obtain the 4-loop 
/3-function of QCD by combining the BPHZ formalism with the algorithm of || for mass- 
less 3-loop 2-point functions in D := 4 — 2e spacetime dimensions. Yet not until very 
recently was a result obtained in [Ifl], and even there the BPHZ formalism was not fully 



exploited. This serves as an indication of the challenge of organizing renormalization. 



In this paper, we exploit the work in [JTTJ, [12], [13], 14], which shows that the joblist of 
renormalization is encapsulated by the coproduct, A, of the Hopf algebra of decorated 
rooted trees and that counterterms are given by its antipode, S. We show how truly 
simple it is to automate renormalization by recursive definitions of A and S and hence 
to calculate the contribution of any diagram to a renormalized Green function, in any 
situation where analytic methods are adequate to evaluate the regularized bare diagrams 
entailed by A. We emphasize that in this work we use only well-tried analytic methods. 
Our aims are to make transparent the Hopf algebra of renormalization, to program it 
efficiently, to apply the program to cases where analysis is possible at large numbers of 
loops, and to report the first fruit of the process of discovery that is thereby enabled. 

In sect. 2 we review 4 crucial formulae from [p], |i~2"f , which capture the entirety of 
renormalization, and rewrite them as a few lines of code in the symbolic manipulation 
language Reduce whose translation to other languages should present little difficulty 
to readers with other preferences. In sect. 3 we give a representation of the operator 
product expansion, recently found by DK [ 14| in the course of extending Chen's lemma [|l(| 



on iterated integrals, and hence providing a powerful test of the correctness of our code. 
In sect. 4 we specialize to the Hopf algebra of undecorated rooted trees, where analytic 
methods and efficient programming allow us to study its 7,813 Feynman diagrams with 
loop numbers n < 12 in a wide variety of field theories and renormalization schemes. 
In sect. 5 we report our first discovery with this new tool, by showing the remarkable 
simplification that results when one combines diagrams with the weights specified by 



Connes and Moscovici 17], in their study of the subalgebra entailed by the diffeomorphism 



group. Sect. 6 offers conclusions and suggestions for further study. 
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Figure 1: The divergences in Tq are generated by 4 skeleton diagrams. 

2 Renormalization by Hopf algebra 



The divergence structure of a Feynman diagram is naturally represented by a tree fl2 



whose vertices (or nodes) represent primitive divergences, with edges connecting vertices 
in a manner that encodes the nesting (or forest structure) of subdivergences. Each prim- 
itive divergence is associated with a skeleton diagram [0, free of subdivergences. To 
distinguish these we need labels, 7*.. Fig. 1 gives an example from QED, entailing one- 
loop skeleton diagrams, 71 and 73, in the electron and photon propagators, and one- and 
two- loop skeleton diagrams, 70 and 72, in the coupling. The assignment of labels is arbi- 
trary, so long as it permits no confusion. Now consider the 6-loop diagram of Fig. 1. How 
shall we encode its divergence structure? 



One obvious method is to form a listQ: 



r 6 = {0,{3},{0,{!},{2}}} 



The first element of the list is a label, which tells us in this case that the subdivergences 
reside inside the one-loop skeleton for the coupling. The rest of the list consists of sublists. 
The first of these is {3}, which encodes the one-loop primitive divergence in the photon 
propagator; the second is {0, {1}, {2}}, which we parse in exactly the same manner as 
before: it tells us that inside a one-loop coupling skeleton there reside a one-loop electron- 
propagator skeleton and a two-loop coupling skeleton. It is also rather convenient that 
this iterated structure, list=head+sublists, allows us to preserve the order of sublists, 
which may later be useful in taking traces over spins. 

1 In the equivalent information was encoded by parentheses, written in the reverse order. For 
present purposes, it is more efficient to work with lists, with the outermost label at the head. 
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Figure 2: Divergences encoded by a rooted tree, corresponding to the list (1). 

An equally obvious, and entirely equivalent, coding is provided by the rooted tree of 
Fig. 2. It grows (downwards!) from the root 0, which has a pair of branches. Each of 
these branches is itself a rooted tree: the first has root 3, but no branches; the second 
has root and branches which are the branchless trees with roots 1 and 2. Clearly, 
there is a trivial translation between the structure tree=root+branches and the structure 
list=head+sublists in ([J). In each case there is a feature which encapsulates the forest 
structure of divergences: each branch is itself a tree; each sublist is itself a list. 

We shall present the structure of the Hopf algebra in terms of trees, roots, and branches 
(which are trees); we shall encode it in terms of lists, heads, and sublists (which are lists). 
In |12|, |13| it is shown that overlapping divergences are accommodated by this iterative 
structure, since they require, at most, a sum of lists of decorations. 



2.1 Coproduct 

Let a G A be a tree, where A is the algebra of trees, with a product that corresponds 
merely to the commutative product of the associated Feynman diagrams. A coproduct 
A(a) G A® A, which is not cocommutative, is defined recursively by |TT], [12] 



A(a)=a®e + id®B+(A(B-(a))) (2) 

where e is the empty tree (evaluating to unity); id is the identity map in A; r is the root 
of a; B~(a) G A is the product of the branches of a; and is the operator that maps 
any element c G A to the tree with root r and branches c. The recursion terminates with 

A(e) = e<g>e (3) 

and is effected by applying the product rule 

A(6) = UkHh) (4) 

to the product of branches b = Y\ k bk obtained by removing the root r of a, and then 
recombining the products on the right of the tensor product with the original root r. 

It was shown in [llj that A(a) generates the joblist of the practical field theorist 
engaged upon the renormalization of a Feynman diagram with divergence structure given 
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by the tree a. On the left of the tensor product reside (products of) subdiagrams requiring 
counterterms; on the right reside diagrams in which these subdiagrams shrink to points. 
Most of us will testify, if pressed, that it is perilously easy to forget where one has reached 
when subtracting sub divergences by hand. On the other hand, a correctly programmed 
recursive procedure, in a suitable symbolic manipulation language, will not lose track. An 
implementation in Reduce is provided by 

for all a,b such that first a=af let A2 a*A2 b=A2 append (a, {b}) ; 
procedure D a; 

Al a+sub(af =f irst a, A2{af }*f or k:=2: length a product D part(a,k)); 



where Al and A2 are operators that hold lists on the left and right of A ® A and af is a 
dummy argument that is replaced by first a after each recursion. The let statement in 
the first line sweeps up products on the right of A® A. In the body of the procedure, the 
list representing the tree is beheaded by starting at k:=2; the product of coproducts of 
sublists is taken; then product terms on the right are reheaded by the original root. This 
generates the joblist of renormalization, which finally reads, in Sweedler notation 

A(a) = E fc 4 1) ®4 2) (5) 
with products of trees on the left, and single trees on the right, of the tensor product. 



2.2 Antipode 

An antipode, S, upgrades the bi-algebra defined by A to a Hopf algebra. For any tree 
a G A, the antipode, S(a) G A, is given, again recursively, by [jO], 

S(a) = -a- E^(4 1} )ai 2) (6) 

where J2'k omits the terms a® e and e ® a of the full Sweedler sum (|T). The encoding is 



procedure S a; sub(Al=S,A2=Al,Al a-D a); 



2.3 Counterterm 



Renormalization consists in applying an operator R, on the left, at each recursion in the 
computation of the antipode. In the momentum (MOM) scheme, -Rmom is merely the 
instruction to replace the external momenta by a set of fiducial momenta. In the minimal 
subtraction (MS) scheme, Rms nullifies external momenta and internal masses and takes 
only the singular terms of the Laurent expansion in e. In either case, we have 

R(R(x)) = R(x) (7) 
R{R{x)R{y)) = R{x)R{y) (8) 
R(x + y) = R(x)+R(y) (9) 
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The counterterm, associated by scheme R to a Feynman diagram with divergence structure 
given by tree a, is [[□], [12 

S R (a) = -R(a) - Tf k R(S R {a^)a { p) (10) 

and is hence encoded by 

procedure S_R a; R(sub(Al=S_R, A2=A1 , Al a-D a)); 

2.4 Renormalized Green function 

The Green-function contribution of a diagram with divergence structure a is 



T R (a) = hmE k S R (a ( k 1) )a [ k 2) (11) 



and is hence encoded by 

procedure G_R a; sub(Al=S_R, A2=A1 ,D a); 

which generates 2 n terms in the case of a tree with n vertices, each representing a primitive 
divergence. The limit e — > is not specified in the procedure, but is guaranteed to exist 
when one provides concrete expressions for the dimensionally regularized bare diagrams, 
as we shall do in sect. 4. 

2.5 Program handsgam . red 

The Reduce program handsgam . red, available^ via FTP, adds 5 features. 

1. A noncom instruction prevents reordering of sublists. This is useful in cases where 
one needs to take account of spin structure. 

2. The linearity of R is declared by let statements. 

3. The loop number of a tree is computed by adding the loop numbers of its primitive 
decorations. This is included as an argument of R, to facilitate Laurent expansion 
in the MS scheme. 

4. Results for A, S R and T R are stored, to speed up batch processing of diagrams with 
sub divergences in common, and tests of the operator product expansion. 

5. A test is made for the arbitrarily decorated tree {1, {2}, {3, {4}, {5}}, {6}}, to ensure 
that the program correctly generates the 64 terms of the BPHZ formalism. 

The reader is invited to generate the 64 BPHZ terms by hand, and then to estimate how 
long it might take to generate 24,213,878 BPHZ terms for the 7,813 diagrams of sect. 4, 
without the assistance of Hopf algebra. 



See ftp : //physics . open . ac . uk/pub/physics/dbroadhu/ha/hareadme . txt for details 
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3 Chen's lemma and the OPE 



Suppose we renormalize in the MOM scheme at some fiducial momenta labelled generically 
by pi and evaluate Tmom at some different external momenta p2- Let Tj v[ Qj kl (a) denote 



the result for a Feynman diagram with divergence structure a. In |14[ it is shown that 

r (i.2) (n] _ v r (i,o) / (i)n f (o,2) . (2), n9) 

1 MOM\ a ) — 2^k L MOMl a fc ) L MOM\ a k ) 

where the sum is over the Sweedler decomposition of the coproduct of a, the momenta po 
are arbitrary, and it is understood that 

r*(6) = Wnibj) (13) 

for a product of trees b = TTj bj, in any scheme R. 



Identity ([12]) is readily obtained in the toy model of |TT[ , which represents the Hopf al- 
gebra of undecorated trees by iterations of products of integrals. In that particular case, it 
corresponds to an extension of Chen's lemma fl6| for iterated integrals. It further extends 
to arbitrarily decorated trees, in any field theory. In |T^] it is construed as a representa- 
tion of Wilson's operator product expansion |1| (OPE), whose distinctive feature is to 



allow one to express matrix elements as sums of products with an arbitrary intermediate 
scale. Here we use (0) as a stringent test of our code: program hamomope . red verifies 
this MOM-scheme OPE. 

There is a further result, for the MS scheme. Consider a Feynman diagram, with di- 
vergence structure a, contributing to a massless two-point function, with a single external 
momentum^ p. Then the MS-renormalized result will be a polynomial, Fusi®, L^), i n 

:= log(/i 2 /p 2 ), with degree equal to the loop number, where /i is the renormalization 
scale of the MS scheme. The MOM-renormalized result will be a different polynomial, 
rMOM^i^i)) m Li := \og(pl/p 2 ), where pi is the fiducial momentum. However, there is 
a strong connection between the sets of MS and MOM polynomials for the terms in the 
Sweedler decomposition of the coproduct A(a), namely 

Tms^ — Li) = Y^k^Msiai^i L^Tmqm^^, —Li) (14) 

with the MOM-scheme results serving to transform MS-scheme results, at momentum p, 
to a MS-scheme result with momentum p\. 

We shall apply ([14]) as a check of our methods in sect. 4. After checking that the 
calculus passes this test, one may economize by computing the MS results at p 2 = fi 2 , 
since their dependence on external momentum, determined by the renormalization group, 
is generated by the MOM-scheme results, by setting L„ = in (114]) . 



4 Hopf algebra of undecorated rooted trees 

We now specialize to the Hopf algebra of undecorated rooted trees, generated by a single 
primitive divergence. From the point of view of pQFT, this is a drastic step: the renor- 

3 We adopt a metric such that p 2 > in the spacelike region. 
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malized Green functions of the theory, at a given order n, are obtained by computing all 
decorations with total loop number n and then combining them with the weights given 
by Wick contractions of the perturbative expansion; little interest might seem to attach 
to a set of Feynman diagrams with only one-loop primitive propagator sub divergences. 

In fact there are (at least) 6 cogent reasons for studying such diagrams. 



1. Analytic results for such bare diagrams are immediately available, in terms of highly 
structured products and quotients of 

T(l -z) = exp (7* + T, s >iC(s)z s /s) (15) 

where z is an integer multiple of e, 7 is Euler's constant, and £(s) := J2k>o^-/k s 
is the Riemann zeta function. Hence we shall be able to obtain explicit results, of 
an analytically nontrivial character, for 7,813 diagrams with loop numbers n < 12, 
entailing the irreducible zeta values {C(3), C(5), C(7), C(9), C(H)} an d powers of n 2 
from, for example, <^(12) = 69l7r 12 /638512875. At the very least, this provides a 
healthy workout for our Hopf algebra code. More importantly it enables us to look 
for structure in the counterterms, lacking in the singularities of bare diagrams. 

2. The 12 pure rainbow diagrams in our set of 7,813 correspond to trees in which no 
vertex has fertility greater than unity. They have a remarkable feature, observed 
in [^0| and proved to all orders in [21, ^] : their MS counterterms are rational 



polynomials in 1/e, bearing no trace of the derivatives of logT(l — z) at z = 0, 
despite the appearance of odd zeta values and powers of 7r 2 at every step of the 
BPHZ formalism. This potent example of the simplicity of counterterms was the 
stimulus for the investigations of |23], |24], |26], |27|, ^8), |29|, |3(J , where we discovered 
remarkable features of the nexus of knot /number /field theory pif . 

Trees in which only the root has fertility greater than zero come from pure chains 
of one-loop self-energy divergences. In the case of QED these generate the ultra- 
violet Landau singularity of the photon propagator. In the case of heavy-quark 



effective theory (HQET) they produce a renormalon structure |32] that frustrates 
unambiguous Borel resummation of the perturbation series. 

We may probe even deeper into the number-theory content of the complex of all 
nestings of rainbows and chains, constituting the full Hopf algebra of undecorated 
rooted trees, by embedding it inside a two-loop diagram, as was done merely for 
chains in |33|. Then we will encounter irreducible multiple zeta values p4| , |35| , [36| , 

37J of the types C(5,3) := E m >n>o 1 / m5 ^ 3 and C(3,5,3) := E fe>m>n>0 l/^m 5 n 3 , 



associated by pQFT to the unique positive 3-braid 8-crossing knot and the unique 
positive 4-braid 11-crossing knot, in the 7- loop analysis of ||23||. 



5. Next, and most excitingly, this seemingly trivial set of Feynman diagrams was shown 
in |12| to provide a solution of a universal problem in Hochschild cohomology. 



6. Finally, and most intriguingly, there is a unique combination of these diagrams at 
loop number n which represents a grading of the commutative part of the Hopf 
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algebra Tix of the diffeomorphism group (T2J. Specifically, if one adds the dia- 



grams with integer weights derivable from work by Connes and Moscovici [17 



one 



expects to see some simpler structure, the details of which neither they nor we an- 
ticipated. Thanks to the discovery machine developed in this section, we shall have 
an interesting result to report in sect. 5, with implications for the relation between 
noncommutative geometry and quantum field theory. 



4.1 Enumeration of rooted trees 



The numbers, N(n), of rooted trees with n vertices form a sequence beginning [38 

1 , 1 , 2 , 4 , 9 , 20 , 48 , 115 , 286 , 719 , 1842 , 4766 , 12486 , 32973 , 87811 , 235381 , 634847 

which grows exponentially, giving 

N(250) = 517763755754613310897899496398412372256908589980657316271 
041790137801884375338813698141647334732891545098109934676 

trees with 250 vertices. Already one sees a nontrivial feature of the iterated structure: 
tree=root+branches, with every branch being itself a tree. The asymptotic growth 

N(n) = ^c n (l + 0(l/n)) (16) 

entails constants 

b = 0.43992401257102530404090339143454476479808540794011 

98576534935450226354004204764605379862197779782334. . . 

c = 2.95576528565199497471481752412319458837549230466359 

65953504724789059647331395749510866682836765813525. . . 

that do not appear to be algebraic numbers, since neither the PSLQ |5jJ nor the LLL [^ 
algorithm was able to find rational polynomials that fit the first 100 digits and correctly 
predict the next 10. The origin of their (presumed) transcendent ality resides in the 
functional equation 

A(x k ) 

log A(x) - A{x) = log x + — V — (17) 

k>i k 

for the generating function A(x) := J2 n >o N(n)x n , which develops a square-root branch- 
point at x = 1/c, near to which 1 — A(x) ~ y/l — ex. Setting x — 1/c we obtain 

logc=l + E A(C ~ fc) 



i, j k 
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which allows one to determine c, by iterative solution, to an accuracy of about l/iV(250) = 
O(10~ 114 ), after determining {N(n) | n < 250} from (|17|), at small x. Then Taylor 
expansion at x = 1/c yields 

2tt6 2 = 1 + £ ^p- (19) 

k>l ° 

with the regularity of the r.h.s. of (O) forcing nonlinear relations between the coefficients, 
c fc , of = 1 - Efe> Cfe(l - cxf^. 

We hope that these asymptotic results may help colleagues to identify the Lie alge- 
bra dual to the Hopf algebra T~Lr of rooted trees, which solves a universal problem in 
Hochschild cohomology [12]. 



4.2 A variety of models 

We shall consider 10 models that generate values for the bare diagrams: 6 of these corre- 
spond to well-defined field theories; the remaining 4 are instructive toy models that probe 
the analytic structure of bare and renormalized diagrams, while ignoring fine details pro- 
duced by spin and gauge dependence. 



4.2.1 BPHZ model 



In [l 1] a BPHZ toy model was developed, by iteration of products of integrals. To each 
tree, a, it assigns a dimensionally regularized bare value, B(a), obtained by the recursion 

B(a) = B n (e)Y[ k B(b k ) (20) 

where n is the number of vertices of a; h k are the trees formed from its branches; and 

BM = (21) 
ne 

where L(e, 5) is regular at both e = and 5 = 0. This structure is generic; it covers all 
the field theories that we study. The BPHZ model is a toy only because of the simplicity 
of its master function, which is 

L BPRZ (e,5) = -^-\- £ (22) 

Sin 7T0 

resulting from nullifying all dimensionful parameters except for an infra-red regulator, A, 
in the outermost integration. 



4.2.2 Heavy-quark model 

More realistically, consider the propagator of a heavy quark in HQET, at virtuality uo < 0. 
Neglecting a rational function of e and 5, which we later restore, the master function is 



r(i- £ )r(i + 2S) 

^ hqM) - r(l-2e + 25) ( -^' 
whose T functions are immediately apparent in |l I 
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4.2.3 Covariant QFT model 



For a massless- quark propagator in QCD, with spacelike momentum p, the T functions 
in give 

Lqft(£ ' 6) ~ T(l-e-6)T(l-e + 6) iP ] (24) 
with a rational factor which will be supplied later. 



4.2.4 MZV model 

To get multiple zeta values |27j (MZVs) in the bare diagrams, we may embed the full 
complex of chains and rainbows in a two-loop diagram, as was done for pure QED chains 
in [|33|| . Here, we embed it in a finite diagram, so as to remain within the undecorated 
Hopf algebra. This corresponds to using the recursive process (P0|, pT] , p^ ) , to compute an 
n-loop bare value for a covariant propagator diagram, and then multiplying by 



where eq. (17) of our paper with John Gracey |26[] gives the general result for the 2-loop 



2-point function S(a, b, c, d), which generates MZVs via a 3 F 2 hypergeometric series. With 
ni+n 2 = Ti3 + «-4, we have S {riis , n 2 e , n^e , n^e) = 3(nin2 — n3n 4 )e 2 £(3) +0(e 3 ), and hence 
obtain irreducible MZVs at level n + 2 in the MS counterterms from undecorated trees 
with n > 6 loops. At n = 9 we probe as deeply into the relation between knots and 
numbers as we did in |[23|| ; at n — 10, we reach the 12-crossing knots of [^, 27j , which 



revealed an unexpected connection between MZVs and alternating Euler sums [35 . 



4.2.5 Field theories 

As already indicated, we have only to supply a rational function of e and S to convert (|23|) 
to the realistic case of a heavy-quark propagator in HQET. In an arbitrary covariant gauge, 
with a gluon propagator g^/q 2 + (x — ^Q^Qu/l 4 , we obtain 

£hqet(£, 5) = 3 ~ 2£ ~_ (1 2 ~ 2£)X X HQ (s, 5) (26) 

which is divergence-free in Yennie gauge, with x = 3. We may compute in an arbitrary 
gauge, and also specialize to 

WMH^^hqM) (27) 
in Feynman gauge, with x — 1> an d 

Lh Q l(£,5) = ^^L HQ (£,5) (28) 
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in Landau gauge, with x — 0. 

Similarly, we can convert (B3j) to the realistic case of a light-quark propagator in 
QCD. The master function vanishes, identically, in Landau gauge. Hence we compute in 
Feynman gauge, with 

w«. o = - (1 _ 2 i 1 :; ) ) ( ( 2 1 :^ s) w«, *> w 

For the fermion propagator in Yukawa theory we have 

^Yuk(£, 5) = ~ ^ _ - ^4 ^ - _ L QFT (e, 5) (30) 

and for 3 theory, in 6 — 2e dimensions, we obtain 

6(1 s\ 

Lphl(£ ' 5) = -(l-e-6)(2-e-5)(3-e-5) LQFT{£ > 6) (31) 
after removing irrelevant multiples of the couplings. 

Thus one chooses a theory, or toy, by setting the model switch to one of the 10 values 
{bphz, hq, qf t, mzv, hqet, hqf , hql, qcd, yuk, phi}, corresponding to the processes ([22}-|3l"|) . 



4.3 A variety of schemes 

We shall investigate the 10 models in 9 renormalization schemes. This multiplicity of 
schemes results from a pair of ternary switches, mscheme and gscheme, each of which may 
take the values 0, 1, or 2. 



4.3.1 Momentum, minimal and nonminimal schemes 

Our first ternary switch, mscheme, chooses between the MOM scheme, the MS scheme, 
and a nonminimal scheme (NMS). 

With mscheme :=0, we retain the full e dependence of the bare diagrams and generate 
counterterms by replacing the external momentum by a fiducial momentum. 

With mscheme : =1, we retain in the counterterms only the singular terms of the Laurent 
expansion, after nullifying dimensionful parameters. 

With mscheme: =2, we retain the finite term as e — > 0, in addition to the singular 
MS terms. Such finite counterterms occur in the delicate handling of 75 in dimensional 



regularization [32, 42]. Like the MOM and MS schemes, this NMS scheme satisfies 



4.3.2 G schemes 

In j|3| it was observed that the appearance of n 2 (but not of 7r 4 ) may be suppressed, 



in covariant field theory, by absorbing suitable V functions into the D-dimensional cou- 
pling, which already absorbs the universal factor l/(Aix) D l 2 . Our second ternary switch, 
gscheme, reflects this freedom. 
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With gscheme:=0, we leave the master function fl2"2] ) as it stands and modify ( |2"B| , |2"4l ) 
only by including the canonical MS factor exp(7e), which suppresses Euler's constant, 7. 

With gscheme:=l, we divide the master functions (]22]-p4D by their values at S = e 
and unit scale, thus obtaining Bi(e) = 1 for the one- loop diagram, when the scale — i.e. 
A in fl22|); — 2u in fl23|); or p 2 in ( |24"D — is set to unity. This suppresses 7 in the HQ and 
QFT cases, and also ((2) = ir 2 /6 in the QFT case. We know from that it 2 is intrinsic 
to HQET counterterms. 



With gscheme:=2, we leave the master function (|22|) as it stands and multiply ( |23| . [24|) 
by T(l — e). This likewise suppresses 7 in the HQ and QFT cases, and also C(2) = n 2 /6 
in the QFT case, yet it leaves a nontrivial value for the one-loop term Bi(e) and hence 
seems less contrived than the G scheme in 



We believe that the 9 schemes entailed by our pair of ternary switches subsume much 
of current practice in pQFT. We make no attempt to implement dimensional reduction, 
analytic regularization, or differential renormalization, though these appear to present no 



problem of principle to the global description of [11, 12 



4.4 Computational strategy 

There are £n=i N(n) = 7, 813 undecorated diagrams with n < 12 loops. To obtain their 
renormalized values we must compute £* 2 =1 2 n N(n) = 24, 213, 878 BPHZ terms. To do 
this in 10 models and 9 renormalization schemes would require us to process more than 
2 x 10 9 BPHZ terms, with nth order polynomials of a logarithm appearing at n loops. 
These polynomials entail 90£* 2 =1 2 n (n + l)jV(n) = 27,804,356,640 coefficients. These 
coefficients contain, in general, products of odd zeta values and powers of it 2 . The number 
of rational coefficients of such products is 90 £n=i 2 n N(n) ££ =0 C(k) = 335, 708, 683, 560, 
where the multiplicity, C(k), of zeta products with levels up to k forms the sequence 



1,1,2,3,4,6,8,11,14,19,24,31,39,49,61,76 



for < k < 15, neglecting the proliferation of the irreducible MZVs of model (|25j), 
which increases these integers to give partial sums of the Padovan sequence |27] . Many of 
these 3.3 x 10 11 rationals will be ratios of integers with O(10) decimal digits. From this, 
it is apparent that a complete analysis, up to 12 loops, might entail processing several 
terabytes^] of exact integer data. Thus computational efficiency is at a premium. 

Clearly, the best strategy for any computation that involves all 7,813 diagrams is an 
Aufbau, in which results at n loops are held in core memory, and then used at n + 1 loops, 
with only a single new recursion of the algorithms (QJ^,^) for the coproduct, antipode 
and bare diagrams. To achieve a 12-loop renormalized result, one needs to hold 13 terms 
from each Laurent expansion, at every loop number n < 12, but need not hold the results 
for the majority of diagrams, namely those 4,766 with n = 12, if only weighted totals, 
such as those presented in sect. 5, are required as output. 

For loop numbers n < 15, we would be talking of petabytes. 
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4.5 Results and timings 



As benchmarks, we present timings for obtaining the renormalized values for all diagrams 
up to a given loop number. As interesting weights, for the output of results, we choose 
those suggested by the work of Connes and Moscovici [|17[] (CM) and derived in sect. 5. 
We emphasize that these weighted results were obtained by summing over the results in 
the full Hopf algebra of undecorated trees, not by specializing to the CM case ab initio. 

Allocating 128MB of core memory to Reduce on a 533MHz DecAlpha machine we 
computed the case model : =bphz ; mscheme : =0 ; gscheme : =0. The 10-loop CM-weighted 
counterterm, at scale A = 1, has finite part 

5 CM,finite = 2140469117T 10 

10 ~ 2112 1 ' 

obtained in less than 2 minutes. The 12-loop result 

5 CM,finite = 435567078937017T 12 
12 "~ 1048320 1 ' 

took less than 2 hours. We would be interested to learn of any method that can obtain 
this result, from explicit computation of all 7,813 counterterms at loop numbers n < 12, 
in a significantly shorter time. 



The next step was to test the code against the MOM-scheme OPE (|T2|), which was 
verified, as were the MS-scheme OPE (TI4T) and its NMS-scheme version 



rNMs( a ) — Li) — X^rNMs( a i \ L^)T M o M (a[, , —Li) (34) 

Thus we may economize, by computing only MOM-scheme results at arbitrary momen- 
tum, while p 2 = fi 2 suffices in the MS and NMS cases, thanks to the renormalization 
group. This helps to keep the memory requirements within practical bounds, since the 
MOM-scheme renormalized Green-function contributions are clearly rational polynomials 
of a log, and are independent of the gscheme switch. 

The modifications of spin and gauge-dependence, in (|26|-pl|), introduce no new analytic 
feature; they merely proliferate terms by mixing products of zeta values with different 
levels. Hence we were content to compute theories (|27H 3T1) to 10 loops. We omitted the 



arbitrary-gauge case fl26|), which generates polynomials of degree n in x at n loops. We 
postpone presentation of results for the MZV model (|25D to a later publication, which 
will address the associated knot theory [EHl El, 12(1 



At 11 loops, we encountered no difficulty in processing models (j22|-p4|) ; at 12 loops, 
128MB of core memory was sufficient only to process model fl2"2"|) . In Table 1 we present 
the times (to the nearest minute) taken to compute all the counterterms and renor- 
malized values up to the specified loop number, at p 2 = fi 2 in the MS scheme, with 
gscheme=l, in 8 models. All timings refer to Reduce3.5, which was allocated 128MB 
of core memory on a 533MHz DecAlpha; changes of scheme have little effect on them. 
Results for the CM-weighted counterterms and renormalized Green functions are in the 
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12 files ha<model><loops> . out. In haphilO.out one encounters, inter alia, a 25-digit 
integer with a 22-digit prime factor. 

Table 1: Time taken to compute all MS counterterms and renormalized values. 



model 


bphz 


hq 


qft 


hqf 


hql 


qcd 


yuk 


phi 


bphz 


hq 


qft 


bphz 


loops 


10 


10 


10 


10 


10 


10 


10 


10 


11 


11 


11 


12 


minutes 


3 


7 


5 


15 


15 


11 


10 


14 


17 


39 


31 


124 



5 Connes-Moscovici weights 



In [|T|] an important connection was found between the Hopf algebra Tix of the dif- 
feomorphism group, studied by Connes and Moscovici |T7| , and the Hopf algebra Hr 
of undecorated rooted trees |IT|, for which we now have 10 representations, generated 



by p2^3T]). The commutative part of H,t is a subalgebra of "Hr. Along with the obvi- 



ous rainbow subalgebra of ||20|| , it exhausts Jl2[ the proper Hopf subalgebras of Hr. In 
consequence, there exist nonzero CM weights for our 7,813 diagrams, such that summing 
rz-loop diagrams with these weights, we arrive at representations of the grading of the CM 
subalgebra, with results for the counterterms that we expect to be as distinctive as for 



the rainbow diagrams in [20 



5.1 Computation of CM weights 

The CM weights are defined as follows. Suppose that one has the CM combination at 
loop number n. To generate the CM combination at loop number n + 1, one grows each 
previous tree to form n new trees, by attaching a new vertex to each of the original 
vertices, in turn. The starting point, at n = 1, is the unique branchless rooted tree, with 
CM weight W({1}) = 1. It is not difficult to encode this algorithm, to obtain the CM 
weights for the 7,813 trees with n < 12, though our first encoding ran for more than 
an hour, until it had finally distributed 11! = 39,916,800 terms between the 4,766 CM 
weights of 12- loop diagrams. The results are in the file hagenhaw. out, with a coding of 
diagrams given in hagendum . out . 

The defining algorithm for these CM weights appears to be very different from those 
for the coproduct and antipode of Hr. For the latter, we could write simple algorithms 
that were blind, at any given recursion, to all structure save that immediately below the 
root (or head) of the tree. By contrast, CM weights appear to feel all the way down to 
the last vertices of the trees, which we refer to as feetf]. An efficient head-first algorithm 
for the CM weight, W(t), was not easy to find. We finally achieved it, as follows. 

The CM definition translates to the feet-first recursion 

W(t)U(t) = V ; in/ ; )lIi/ ; ) (35) 

5 Recall that a mathematical tree has, like a family tree, its root at the top. 
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Figure 3: A 12-loop diagram based on a one-loop skeleton. 




Figure 4: The divergence structure of Fig. 3, exhibiting permutation symmetry. 

where the fj are all the trees obtained by removing a single foot from t, and II(t) is 
a permutation factor, equal to the number of trees that are indistinguishable from t by 
permutations of branches originating from any vertex. For example, the 12-loop Feynman 
diagram of Fig. 3 has a permutation^ factor n(t 12 ) = 2! x 3! x 3! = 72, as is apparent 
from its divergence structure 

t 12 := {1, {1, {1}, {1}, {1}}, {1, {1}, {1}, {1}}, {1, {1, {1}}}} (36) 

illustrated in Fig. 4. 

Now the remarkable thing about the defining formula (|35| ) is that it produces no 
prime factor greater than n — 1, for a tree with n vertices. This is not at all clear from 
the formula, since the sum seems to allow the possibility of adding, let us say, the integers 
9 and 4, to obtain the forbidden prime 13. To show that this cannot happen, we define 
another construct: the tree factorial, t\ Its recursive definition is 

t l = nU k b[ (37) 

where n is the number of vertices of t and b k are its branches. In the pure rainbow case, 
where no vertex has fertility greater than unity, t l is the ordinary factorial of the number 
of vertices, n. In the case of Fig. 4, we have t 12 = 12 X 4 X 4 x 3! = 1152. 

6 This permutation factor is quite distinct from any field theoretic symmetry factor. 
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So now we have a feet-first sum for the CM weights (|3"5j), and a head-first product 
for the tree factorial (137T). The trick is to convert the latter to a feet-first sum, and then 



express the former as a simple head-first product. In |L4j it is proven, by induction, that 
the tree factorials b l k of the beheadings of any tree are related to the tree factorials fj of 
its befootings by a wonderfully simple formula: 

n^=Ei (38) 

k °k j Jj 

which produces, for example, the equality 

111 6 1 

_ x _ x = I (3Q) 

4 4 3! 11x4x3x3! 11x4x4x2 v; 
in the case of Fig. 4. Then a second induction gives 

Wit) = ^_ (40) 

with t ! given by (|37D, and TI(t) likewise given by a head-first recursion, namely 

U(t) = n(t)Umk) (41) 
k 

where ir(t) is the number of indistinguishable permutations of the branches of t. 

This reduces the computation of the CM weight W(t) to a rather simple procedure. 
One associates a pair of integers to each vertex of a tree: the first is the number of vertices 
of the subtree with this root; the second is the number of indistinguishable permutations 
of its branches. For a tree with n vertices, one obtains W(t) from (f40|) by dividing n\ by 
all these integers. In the case of Fig. 4, 

12! 

W(t 12 ) = — — — - — - - = 3x5 2 x7xll = 5775 (42) 

v ; (12 x 4 x 4 x 3!) x (2! x 3! x 3!) K J 

gives the rather nontrivial Connes-Moscovici weight of the Feynman diagram of Fig. 3. 
Computing the CM weights of all 4,766 12-loop diagrams, and combining results from 
sect. 4 with these weights, we expect to find something as interesting as the rational 



rainbow results of 2C . 



5.2 Results with CM weights 

We do indeed find rationality. The Laurent expansion of the 12-loop CM antipode of the 
BPHZ model (E3), at A = 1, is 



. rx , 155925 43556707893701tt 12 

s ™ = msm — + 0(£ » (43) 
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with no singular term that entails tt. As far as Connes and Moscovici are concerned, the 
Bernoulli numbers generate no subleading Laurent residues! Moreover this depends upon 
no property of the Bernoulli numbers. With arbitrary coefficients, Ck, in the expansion 



B n (e) = -ll + Y / C k (nne) 2k ) (44) 

ne V k>0 J 

the 12-loop CM antipode would have the same 7r-free singularity. Yet the finite part of 

Sg M = 155925 - |51321600C 6 + 21837600C 5 Ci + 7001856(^2 + 3232224C 4 C? 

+ 2534490C 3 2 + 2208576C 3 C 2 Ci + 229020C 3 C , 1 3 + 158336C 3 

+ 146550C?C? + 9180tf 2 C? + ^C^tt 12 ) + 0(e 2 ) (45) 

shows the nontrivial processing of coefficients that results from the recursions for the 
coproduct (H), the antipode (H), and the bare diagrams (^Oj), in the presence of the CM 
weights (gD|). 

With CM weights, the sole singular n-loop counterterm of the BPHZ model ( p2|) is of 
a rational, diagonal, quadratic form, W2(n)/(— e) n , where 

w 2 (n) ■.= £ ^fl = i £ w(t fc )n(t fc )w(**) = ^Vr 1 ( 46 ) 
k=l l k n - k=l z 

posits a closed form, for the sum over all trees tk with n vertices, which we have verified 
for n < 12, yet have not proved in general. Moreover, we found that W / 2(n)/(— e) n is also 
the sole n-loop MS counterterm in the CM- weighted QFT model of (^) at n < 11 loops. 
The minimal subtraction scheme is such common practice in pQFT because it retains in 
counterterms only those (products of) zeta values that are strictly necessary to render 
the renormalized Green functions finite. CM weights combine Feynman diagrams so as to 
annihilate all zeta values in model (p4|), just as they annihilate powers of tc 2 in model (|22|). 

Inspection of (^) and (^4j) shows that each is an even function of 5. Thus we investi- 
gated the most general Ansatz with this property, namely 

B n (e) = — E C hk {ne) 2 ^e k (47) 

n£ j,k>0 

The MS counterterm from an individual 10-loop diagram involves up to 302 terms, formed 
from products of {Cj y k \ < 2j + k < 10}. It took an hour to compute all 719 10-loop 
counterterms, with symbolic values for the coefficients in (|47|). Combining them with CM 
weights, we annihilated 301 products of the arbitrary coefficients, and were left with only 
the expected multiple, W 2 (10) = 9!/2 9 = 2835/4, of (-C 0fi /e) 10 . 

We find it significant that the ineluctable divergences of renormalization are so re- 
sponsive to ideas from noncommutative geometry. This strengthens our opinion that 
renormalization is of real interest, in its own right, and far from being a cause for regret. 
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6 Conclusions and prospects 



In this paper, we automated the process of renormalization. Its Hopf algebra structure, 
imposed by the necessity to generate local counterterms, was summarized in a few lines of 
code. The generality of this algebraic structure extends to the renormalization of realistic 
particle physics problems, as in the Standard Model. Here, we computed deep into the 
perturbation expansion in cases where the remaining analysis was straightforward. 

Different renormalization schemes and changes of scale were implemented with ease, 
thanks to the convolutions ( |i2| , P] ). We regard no renormalization scheme as better defined 



than others; all can be treated on the same algebraic footing. Indeed, it will be shown 
elsewhere [TJ], WM that a minimal subtraction scheme can, for example, be obtained as a 



BPHZ scheme, whose renormalization scales are indexed by rooted trees. 

The results in sect. 5 indicate a deep connection between QFT and noncommutative 
geometry. There is an index-theoretic flavour, in the annihilation of zeta values by Connes- 
Moscovici weights, which deserves further study. In particular, the new relation (f40|), 
between tree factorials and CM weights, will be of great help in making the connection 



between counterterms and diffeomorphisms, as was envisaged in |L2| and will be made 
more precise in f45|. We also hope that the asymptotic enumeration of sect. 4.1 will help 
Lie-algebra specialists to solve the outstanding problem of identifying the dual of Hr. 

Our construction of a renormalization engine demonstrates that the remaining chal- 
lenges in computational pQFT lie in two related areas: elucidation of algebraic rela- 



tions ||!| f|7| between decorations; analysis of those iterated integrals jnj, or multiple 
sums [H, that are the concrete representations of the truly primitive elements which sur- 
vive this filtration. 

We envisage that both endeavours will be illuminated by an interplay [|1^] with non- 



commutative geometry |17 



Acknowledgments: We thank the ESI for generous hospitality in Vienna; Alain Connes, 
Victor Kac and Ivan Todorov for discussions there; Bob Delbourgo and John Gracey for 
collaborations that contributed to this work; and HUCAM for grant CHRX-CT94-0579, 
which contributed materially to our meeting over the last 4 years. 



18 



References 

N.N. Bogoliubov, 0. Parasiuk, Acta Math. 97 (1957) 227. 

K. Hepp, Comm.Math.Phys. 2 (1966) 301. 

W. Zimmermann, Comm.Math.Phys. 15 (1969) 208. 

J. Collins, Renormalization (Cambridge Univ. Press, 1986). 

G. 't Hooft, M. Veltman, Nucl.Phys. B44 (1972) 189. 

K.G. Chetyrkin, F.V. Tkachov, Nucl.Phys. B192 (1981) 159; 
F.V. Tkachov, Phys.Lett. B100 (1981) 65. 

D.J. Broadhurst, Z.Phys. C54 (1992) 599. 

L.V. Avdeev, Comput.Phys.Commun. 98 (1996) 15. 

D.J. Broadhurst, |hep-th/980309l| , to appear in Eur.Phys.J. C; 
[hep-th/ 9805025 , to appear in Eur.Phys.J. C; 



|hep-th/9806T7l , submitted to Eur.Phys.J. C. 



T. van Ritbergen, J.A.M. Vermaseren, S.A. Larin, Phys.Lett. B400 (1997) 379, 
|hep-ph/ 9701390 . 



D. Kreimer, Adv.Theor.Math.Phys. 2 (1998) 303, |q-alg/970702 ( o . 



A. Connes, D. Kreimer, Hopf Algebras, Renormalization and Noncommutative 
Geometry, |hep-th/9808042| , to appear in Comm.Math.Phys. 

D. Kreimer, On overlapping Divergences, |hep-th/9810022| . 

D. Kreimer Chen's iterated Integral represents the Operator Product Expansion, in 
preparation. 

A.C. Hearn, REDUCE User's Manual, Version 3.6, RAND publication CP78, July 
1995. 

T.K. Chen, Bull. Amer. Math. Soc. 83 (1977) 831. 

A. Connes, H. Moscovici, Hopf Algebras, cyclic Cohomology and the transverse 
Index Theorem, IHES/M/98/37, |math.DG/9806109| . 

J.D. Bjorken, S.D. Drell, Relativistic Quantum Fields (McGraw-Hill, New York, 
1965). 

K.G. Wilson, Phys.Rev. 179 (1969) 1699. 



D. Kreimer, Phys.Lett. B354 (1995) 117, [hep-th/9503059 



19 



R. Delbourgo, A. Kalloniatis, G. Thompson, Phys.Rev. D54 (1996) 5373, 
|hep-th/ 960510^ . 

R. Delbourgo, D. Elliott, D.S. McAnally, Phys.Rev. D55 (1997) 5230, 
|hep-th/961lT50| . 



D.J. Broadhurst, D. Kreimer, Int. J.Mod.Phys. C6 (1995) 519, |iep-ph/9504352 



D.J. Broadhurst, R. Delbourgo, D. Kreimer, Phys.Lett. B366 (1996) 421, 
|hep-ph/9509296 . 

D. Kreimer, J.Knot Th.Ram. 6 (1997) 479, |q-alg/9607022j 

D.J. Broadhurst, J.A. Gracey, D. Kreimer, Z.Phys. C75 (1997) 559, 
|hep-th/9607T74 . 



D.J. Broadhurst, D. Kreimer, Phys.Lett. B393 (1997) 403, |hep-th/9609T28 
D. Kreimer, Eur.Phys.J. C2 (1998) 757, |hep-th/9610T28| . 
D.J. Broadhurst, A.V. Kotikov, |hep-th/9612"013l to appear in Phys.Lett. B. 
D. Kreimer. Acta Phys.Polon. B29 (1998) 2865, |hep-th/9807i25| . 
D. Kreimer, Knots and Feynman Diagrams (Cambridge Univ. Press, in press). 
D.J. Broadhurst, A.G. Grozin, Phys.Rev. D52 (1995) 4082, |hep-ph/9410240| . 
D.J. Broadhurst, Z.Phys. C58 (1993) 339. 

J.M. Borwein, D.M. Bradley, D.J. Broadhurst, Elec. J. Combinatorics. 4 (1997) R5, 
|hep-th/ 961 lOOl . 



D.J. Broadhurst, [hep-th/96041~28| , to appear in J.Math.Phys. 

J.M. Borwein, D.M. Bradley, D.J. Broadhurst, P. Lisonek, Special Values of 
Multidimensional Poly logarithms, CECM report 98-106, May 1998, submitted to 
Trans.A.M.S. 

J.M. Borwein, D.M. Bradley, D.J. Broadhurst, P. Lisonek, Combinatorial Aspects 
of Euler Sums, CECM report 98-107, June 1998, to appear in Elec. J. Combinatorics. 

N.J. A. Sloane, On-line Encyclopedia of Integer Sequences, sequence A000081, 
http : //www. research, att . com/p hi as/sequences; 

F. Bergeron, G. Labelle, P. Leroux, Combinatorial Species and Tree-like Structures 
(Cambridge Univ. Press, 1998) p. 279. 

[39] H.R.P. Ferguson, D.H. Bailey, S. Arno, Analysis of PSLQ, an integer relation 
finding algorithm, NASA-Ames Technical Report, NAS-96-005, to appear in 
Math.Comp. 



20 



[40] A.K. Lenstra, H.W. Lenstra, L. Lovasz Math.Ann. 261 (1982) 515; 
R. Kannan, A.K. Lenstra, L. Lovasz, Math.Comp. 50 (1988) 235. 

[41] D.J. Broadhurst, A.G. Grozin, Phys.Lett. B267 (1991) 105; 
Phys.Lett. B274 (1992) 421. 



[42] D.J. Broadhurst, A.L. Kataev, Phys.Lett. B315 (1993) 179, |hep-ph/ 9308274 



[43] K.G. Chetyrkin, A.L. Kataev, F.V. Tkachov, Nucl.Phys. B174 (1980) 345. 

[44] R. Delbourgo, D. Kreimer, in preparation. 

[45] A. Connes, D. Kreimer, in preparation. 

[46] D. Kreimer, J.Knot Th.Ram. 7 (1998) 61, |hep-th/96120T0 . 



[47] D.J. Broadhurst, D. Kreimer, Phys.Lett. B426 (1998) 339, |iep-th/96120TT 



21 



